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LONG-TERM  GOALS 

The  long-term  goals  of  this  project  are: 

(1)  to  develop  a  theoretical  framework  to  quantify  turbulence  induced  NPZ  interactions. 

(2)  to  apply  the  theory  to  develop  parameterizations  to  be  used  in  realistic  environmental  physical 
biological  coupling  numerical  models. 

APPROACH 

Connect  the  Goodman  and  Robinson  (2007)  statistically  based  pdf  theory  to  Advection  Diffusion 
Reaction  (ADR)  modeling  of  NPZ  interaction. 

Background 

A  nonlinear  model  for  biological  and  physical  dynamical  interactions  in  a  laminar  flow  field  being 
upwelled  into  the  mixed  layer  Robinson(1999)  (Fig  1,  below)  has  been  extended  to  turbulent  flow 
(Goodman  and  Robinson,  2007).  The  approach  of  the  Goodman  and  Robinson  theory  has  been  to 
develop  a  probability  density  function  (pdf)  for  the  turbulent  displacement  field  and  use  that  to 
calculate  the  turbulence  induced  biological  interaction  (TIBI)  terms,  i.e.  <PiNj>,<PiZJ>,<ZiNj>, 

where  N  i,Pi,Zi  are  the  ith  component  of  a  field  of  different  nutrients,  phytoplankton,  and  zooplankton 
embedded  in  the  turbulent  field.  Contrast  the  TIBI  terms  with  the  biological  turbulent  flux  terms 

<  u '  /V,  >,<u'  P  >,<u'Zj  > .  The  formalism  for  modeling  the  latter  type  of  terms  are  well  developed, 
typically  involving  some  type  of  eddy  diffusivity  or  higher  order  closure  such  as  Mellor  and  Yamada, 
(1982).  However,  at  present,  no  biodynamical  basis  for  closure  of  the  TIBI  terms  has  been  developed. 

The  current  approach  to  handling  the  TIBI  terms  is  either:  (1)  to  ignore  them  by  setting 

<  (..),(••),  >=<  (••),  ><  (••),  >  in  an  advection  diffusion  reaction  (ADR)  equation  approach  (Donaghay, 

and  Osborn  ,  1997),  or  (2)  to  perform  a  numerical  simulation  for  the  turbulent  displacement  field  and 
explicitly  calculate  the  TIBI  terms.  The  former,  as  we  will  show  below,  can  typically  result  in  a  large 
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overestimate  of  the  effect  of  turbulence  on  these  interactions,  while  the  latter  is  very  limited  in  its 
domain  size,  suffers  from  difficulty  in  imposing  realistic  boundary  condition  at  the  laminar  turbulent 
interface,  and  can  only  reveal  significant  physics  of  the  TIBI  terms  with  a  large  number  of  repetitive 
runs,  which  prohibit  the  size  of  the  computational  domain  to  very  limited  environmental  scenarios.  In 
addition,  turbulence  numerical  models  such  as  LES  and  DNS  are  difficult  to  embed  in  larger 
regional  scale  biophysical  coupling  models.  What  is  needed  to  be  useful  in  the  larger  scale  physical 
biological  coupling  models  is  development  of  a  realistic  parameterizations  of  the  TIBI  terms, 
analogously  to  the  development  of  turbulent  flux  parameterizations  used  in  regional  and  large  scale 
ocean  circulation  models.  In  addition  development  of  such  parameterizations  will  lead  to  new 
physical/  biological  insights  into  the  role  of  the  TIBI  terms. 

WORK  COMPLETED 

In  the  past  year  and  half,  since  the  beginning  of  this  project,  in  order  to  develop  parameterizations  for 
the  TIBI  terms,  we  have  been  examining  the  relationship  of  the  Goodman  and  Robinson  pdf  theory 
approach  to  that  of  the  ensemble  averaged  ADR  approach.  We  have  made  a  major  breakthrough  and 
have  obtained  an  exact  solution  to  the  ADR  equation  for  realistic  boundary  conditions  using  the 
Goodman  and  Robinson  pdf  theory  approach.  With  these  results  we  are  now  in  a  position  to  formally 
examine  various  parameterizations  of  the  TIBI  terms.  Preliminary  results  indicate  that  for  the  single 
component  bilinear  NP  interaction  problem,  the  potential  suitability  of  a  linear  parameterization. 
Below  we  outline  the  results  to  date. 


Consider  first  the  Advection  Diffusion  (AD)  equation  with  no  biological  interaction.  Let  p.  be  some 
biological  scalar  such  the  total  biomass  density  being  upwelled  into  a  turbulent  optically  inactive  mixed  layer  , 
with  p  =  P  +  N  ,  where  P,  N  are  the  phytoplankton  and  nutrient  mass  densities,  respectively.  Using  the 
simple  linear  upwelling  field  of  Goodman  and  Robinson  (2007)  (Fig.  1 )  and,  by  symmetry,  assuming  that 

/v 

turbulent  mixing  is  dominated  by  vertical  (ID)  fluxes  ,  it  is  straightforward  to  show  that  P  satisfies  the  one 
dimensional  AD  equation 


=  0 


with  the  boundary  conditions 

(2a)  —  =  0  at  z  =  0 

dz 


(2b) 


dp 

ap  +  K — c- 

dz 


=  ap0 


atz  =  1 


where  the  depth  z  is  normalized  by  the  mixed  layer  depth,  D;  a 
rate, - formalized  by  the  nutrient  uptake  time,  T  te  ; 

^ ^ upwelling 


— uptake  js  upwelling  strain 

^upwelling 

kY 

k  =  — ,  is  the  normalized  turbulent 
D 
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diffusivity  ,  k  ,  p0  is  the  biomass  density  being  upwelling  into  the  turbulent  mixed  layer,  which  is 
located  at  z  <1 .  Boundary  condition  (2a)  ensures  that  no  material  is  fluxed  out  of  the  mixed  layer  at 
the  surface,  while  boundary  condition  (2b)  ensures  that  mass  flux  is  conserved  at  the  base  of  the  mixed 
layer,  ft  should  be  noted  that  in  that  the  model  pdf  discussed  in  Goodman  and  Robinson  (2007)  once  a 
particle  entered  the  mixed  layer  it  is  trapped  in  the  mixed  layer.  This  lead  to  a  no  flux  boundary 
dp 

condition,  —  =  0  ,  at  the  base  of  the  mixed  layer  and  resulted  in  a  build  up  of  the  integrated  primary 

dz 

production.  Boundary  condition  (2b)  is  a  more  realistic  characterization  of  the  role  of  turbulence  at  the 
base  of  the  mixed  layer,  i.e.  transition  zone  of  Fig  1 .  Eq.  (1)  with  boundary  conditions  (2a)  and  2b)  can 
be  solved  exactly  where  it  can  be  shown  that 


(3)  P  =  Ad  -  Sill  T„G„,  exp[-(P^Z  ]  exp (~aym  (t  - 1„ )] 

where  ym ,  Gm ,  are  eigenvalues  and  eigenfunctions  solutions  of  the  equation 

1  2  p2 

(4)  Pe[y - Pe—]G  +— 'y~  =  0 

y  >  um  2  4  J 

where  Pe  is  the  Peclet’  number  defined  by 


K 


which  measures  the  relative  importance  of  advection  to  turbulent  diffusion.  The  time  parameter  to  is  the 
initial  time  when  turbulence  is  turned  on.  From  the  orthononnal  property  of  eigenfunctions,  it  is 
straightforward  to  show  that 


L  pP 

A„  =  J  dz ' G,„  (z ')  exp(— ^ z 12 ) 
0  4 

Note  that  at  t  =  t0  f)  —  0  Z<1 
We  can  formally  write  (3)  as 


J  dt'H(t0,t')( 

—00 


with  the  Heaviside  function  defined  by 
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t’>t0 

t'<t0 

/\  - 

Eq  (5)  can  be  interpreted  as  the  prescription  for  obtaining  the  average  density  P  —  P  =<^  P  ^ 
from  the  initial  density  p  =  p0H(t0,t')  using  the  probability  density  function  pdf,  F,  namely 

t 

(6)  p=\dt'p(t')F(t'-t) 

—GO 


with 


(7)  Fit t)  =  —  (- =  «Xi: r„A„G„,  exp[-(jPg)Z  ] exp {-aym (t  - 1 ')] 


Pa  dt 


Changing  variables  to  z=exp(-a(7  -t')  we  find  that 

i 

(8)  p  =  |  dzp(z)F(z) 
o 

In  general  we  find  that  an  average  value  can  be  written  as 

t 

<...>=  f 

—00 

Consider  phytoplankton  of  density,  P,  and  Nutrients  of  density,  N,  being  upwelling  into  a  turbulent 
mixed  layer.  Let  the  interaction  of  P,  N  be  a  simple  bilinear  form 

dP 

(9  a)  —  =  PN 
dt 

dN 

(9b)  ——  =  -PN 
dt 

Using  the  pdf,  F,  we  can  write 
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t  1 

(10)  P  =  f  dt'P(t')F(t')  =  f dt'P{z)F{z ) 

—00  0 

t  1 

(1 1)  <  PN  >=  f  dt  P(t  ')N(t  ')F(t  ')=\dt  P{z)N{z)F{z) 

—oo  0 

Consider  now  the  equation  for  P  ,  from  (9a)  and  (b).  It  follows  that 


r>+JV0exp[-(/-/0)}  H<f,h 

which  upon  using  in  (10)  and  (11),  noting  that  N  =  1-P  ,and  with  some  algebra  yields  the  differential 
equation 


dP  dP  d2P 

(13)  --az- — K—r  =  aP0S(z-l)+<PN> 


dt 


dz  dzz 


subject  to  the  boundary  conditions 


(14  a) 


dP 

dz 


=  0  at  z  =  0 


dP 


(14 b)  aP  +  K —  =  0  at  z  =  1 

dz 


Note  that  the  first  term  on  the  RHS  of  (13)  represents  a  source  term  at  z  =1  and  thus  P  can  be  cast  in 
the  exact  same  form  as  equations  (1)  and  (2)  or 


dP  dP  82P 

(15) - az - k — —=<PN> 


dt 


dz  dz z 


(1 6a) 


dP 

dz 


=  0  at  z  =  0 


dP 


(1 6b)  aP  +  /c —  =  aP  at  z  =  1 

dz  0 
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This  proves  that  using  equation  (7)  in  (10)  satisfies  the  ADR  equation,  (15)  subject  to  boundary 
conditions  (16a)  and  (16b),  which  is  the  same  as  (2a)  and  (2b). 

We  are  now  in  a  position  to  use  this  rigorous  approach  to  compare  various  parameterizations  of  the 
TIBI  terms.  Consider  the  bilinear  NP  interaction  discussed  above.  The  simplest  parameterization  is  to 
take  the  relationship  for  induced  turbulent  TIBI  term  as  <  P'N'  >  =  -CPN  ,  which  results  in  the  RHS 
of (15)  as 

<  PN  >  =(1  -C)PN . 

The  minus  is  chosen  since  for  this  case  our  analysis  has  shown  that  <P' N'  >  is  always  negative. 
Note  that  the  standard  approach  (Donaghay,  and  Osbom,  1997)  is  to  set  C  =0  .  Table  I  below  shows 

i 

vertically  integrated  values  of  C,  namely  <  C  >s=  J  dzC  for  different  values  of  a  with  Pe  =  1,  i.e. 

o 

advection  and  turbulence,  playing  an  equal  role.  Also  shown  are  the  exact  solutions  for  <  C  >5 


obtained  with  our  new  pdf  theory,  described  above.  Although  it  should  be  noted  that  we  are  not 
comparing  the  actual  local  values  of  C,  i.e.  C  =  C(z),  the  closeness  between  theory  and  this 
parameterization  approximation  is  very  encouraging.  Note  as  a  increases  the  uptake  time  increases 
relative  to  the  advection  time  and  turbulence  plays  a  decreasing  role  in  TIBI  interaction.  Conversely 
for  small  a  and  thus  small  uptake  time  relative  to  the  advection  time  the  role  of  turbulence  on  the  TIBI 
tenn  increases. 

As  of  this  writing,  I  am  in  the  process  of  writing  this  up  for  publication  either  to  Proc  Roy  Soc  as  a 
follow  up  of  our  previous  manuscript  or  to  JMS. 
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Figure  1.  Upwelling  flow  field  into  an  optically  and  biologically  active  mixed  layer  used  in  the 
Goodman  and  Robinson  (2007)  biodynamical  model.  Note  as  Diner  eases  or  as  the  uptake  time 

increases  relative  to  the  advection  time 


<P'N'> 
Values  of  C  =  — = — 
PN 


Table  I 

for  Pe  =1  and  difference  values  of  a, 
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